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We study the diffusion of tagged hard core interacting particles under the influence of an external 
force field. Using the Jepsen line we map this many particle problem onto a single particle one. 
We obtain general equations for the distribution and the mean square displacement {{xt)'^) of the 
tagged center particle valid for rather general external force fields and initial conditions. A wide 
range of physical behaviors emerge which are very different than the classical single file sub-diffusion 
{{xtY') ~ t^^'^ found for uniformly distributed particles in an infinite space and in the absence of 
force fields. For symmetric initial conditions and potential fields we find {{xt)"^) = ''^2n'J^^ where 2N 
is the (large) number of particles in the system, 7?. is a single particle reflection coefficient obtained 
from the single particle Green function and initial conditions, and r its derivative. We show that 
this equation is related to the mathematical theory of order statistics and it can be used to flnd 
{{xtY) even when the motion between collision events is not Brownian (e.g. it might be ballistic, 
or anomalous diffusion). As an example we derive the Percus relation for non Gaussian diffusion. 

PACS numbers: 05.60.k, 02.50.Ey, 05.40.Jc, 05.70.Ln 



I. INTRODUCTION 

Systems of particles governed by stochastic dynamics 
and exclusion interactions have been studied for a long 
time [iL'sj . One aspect of this problem is the motion of a 
tagged particle sometimes called the tracer particle 0-Q] ■ 
The diffusion of a tagged particle, in a one dimensional 
system of Brownian particles, interacting via hard core 
interaction is a model for the motion of a single molecule 
in a crowded one dimensional environment such as a bi- 
ological pore or channel [B4iHj and in experimentally 
studied physical systems such as zeolites |12l ] , confined 
colloid particles [l3l - [l5| and charged spheres in circular 
channels [13] ■ Since particles do not pass each other such 
diffusion processes are called single file diffusion. 

Confinement of a tagged Brownian particle, due to 
its interaction with other Brownian particles, leads to 
a slow-down of the diffusion of the tagged particle 
{{xt)^) oc t^^^ instead of normal diffusion ((xt)^) oc t 
when the particles in the system are uniformly dis- 
tributed 0,1^1 ■ Such many body problems can be treated 
using methods which exploit the relation between the dy- 
namics of the interacting tagged particle and the motion 
of a particle free of interactions [1, d, [l7l-l20| . In recent 
years at-least six new directions of research emerged, (i) 
The effect of an external field acting on all the particles 
[m or on the tagged particle only [H, HI] has attracted 
attention since pores induce entropic barriers [13] and 
are generally inhomogeneous. In this category the ex- 
amples of single file motion in a periodic potential [23] 
or in a box [25] were considered in detail, (ii) Initial 
conditions may have a profound impact on diffusion of 
the tagged particle. For example if particles start as a 
narrow Gaussian packet the diffusion of the tagged cen- 
ter particle increases linearly in time ((xt)^) c>c t [IH] 
(see details below). Power law initial conditions induce 



{(xtY) oc and ^ is neither I nor 1/2 [13] (see further 
details in text). When an external field is acting on the 
system, it is important to consider non uniform initial 
conditions where the density of particles is determined 
naturally from Boltzmann's distribution, (iii) In some 
investigations the underlying motion is not normal diffu- 
sion, instead the particles inay be non-Brownian, and fol- 
lowing anomalous kinetics [23, [13]. (iv) Usually hard core 
interactions are considered, although the hard problem 
of more general interactions has been recently treated in 
[13, [29I . [30| . Screened hydrodynamic interactions which 
seem important at short times were investigated in [3l| 
both theoretically and experimentally. Granular single 
file diffusion with inelastic collisions shows the typical 
t^/"^ sub-diffusion [13, H^- (v) Interacting particles in 
systems with quenched disorder is yet another challenge. 
The motion of a tagged particle was recently treated in 
the context of single file diffusion in the Sinai model 
[131 ■ (vi) Finally, if the motion of the tagged particle 
is not normal Brownian motion, what stochastic theory 
replaces the usual Brownian-Langevin framework? In 
this direction interesting connections to fractional cal- 
culus emerged [13, 113, [13] . Roughly speaking and under 
certain conditions half order time derivatives (c?^/^/di^/^) 
enter in the Langevin equation [13] and fractional Brow- 
nian noise replaces white noise. 

Here we provide a very general theory of single file dif- 
fusion of the center tagged particle, valid in the presence 
or the absence of an external potential field V{x), as well 
as for thermal and non-thermal initial conditions. Our 
main results reproduce previously obtained formulas and 
many new ones, by mapping the many particle problem 
onto a single particle model. Our method, explained in 
Sec. ni exploits the theoretical concept of the Jepsen 
line [5|, is limited to hard core point particles, but, as 
we show in Sec. |Vl is not limited to Brownian particles. 
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After providing general results in Sec. IIIII we limit our 
attention to symmetric initial conditions and potentials 
where the tagged particle has no average drift. A general 
relation between the mean square displacement of the 
tagged particle and reflection probability of the non inter- 
acting particle is given in Eq. (|37l) . Detailed calculations 
of the mean square displacement of the tagged particle 
then follow, for special choices of force fields and initial 
conditions, in Sec. |TVl As we discuss in sub-section lIIID[ 
in certain limits our problem is related to order statistics 
[srj . a fact worth mentioning since it allows us to solve 
our problem and related ones using known methods. In 
Sec. |Vl we discuss non-Brownian kinds of motion and 
the Percus relation. A brief report of part of our results 
was recently published [2l| . 



II. MODEL AND METHODS 

In our model, 2A^+1 identical point particles with hard 
core particle-particle interactions are undergoing Brow- 
nian motion in one dimension, so particles cannot pass 
each other. The diffusion constant of particles free of 
interaction is D. As mentioned, towards the end of the 
paper we discuss the more general case where the dynam- 
ics between collision events is not necessarily Brownian. 
An external potential V{x) acts on the particles. The 
system stretches from —L to L; however, unless stated 
otherwise we will let L — >■ cx) and obtain a thermody- 
namic limit where N / L is fixed. We tag the center parti- 
cle, which clearly has N particles to its left and N to 
it right. Initially the tagged particle is at the origin 
a; = 0. The motion of a single particle in the absence of 
interactions with other particles is described by a single 
particle Green function g{x, Xq, t), with the initial condi- 
tions (7(0;, xo, 0) = (5(x — xo). In the case of over damped 
Brownian motion the Green function is the solution of 
the Fokker-Planck equation [ssl ] 



dg{x,xo,t) 
dt 



D 



dx"^ 



1 d 

khT dx 



Fix) 



g{x,xo,t), (1) 



where F{x) — ~V'{x) is the force field, T is the tem- 
perature and fcf, is Boltzmann's constant. The initial 
conditions of N particles residing initially to the right 
(left) of the test particle are drawn from the probability 
density function (PDF) Jr^Xq) (/L(a;o)) respectively. We 
consider an ensemble of trajectories and average over tra- 
jectories and initial conditions (see details below). Our 
goal is to obtain the PDF of the position xt of the tagged 
center particle P{xt) in the large N limit. 



A. The Jepsen Line 

A schematic diagram of the problem is presented in 
Fig. [T] for particles in a box. Initial positions of particles 
are given by Xq where j = —N, • • • , 0, • • • TV where j is 




FIG. 1: Schematic motion of Brownian particles in a box, 
where particles cannot penetrate through each other. The 
center tagged particle label is 0, its trajectory is restricted by 
collisions with neighboring particles. The straight line is the 
Jepsen line, it follows vt as explained in the text. As explained 
in the text, in an equivalent non-interacting picture, we allow 
particles to pass through each other, and at time t we search 
for the position of the particle which has A'^ particles to its 
right and to its left (i.e. the center particle). 



the label of the interacting particles (see Fig. [T]). The 
tagged particle whose coordinate is denoted with XT(i), 
is the center particle j = (bold line in Fig. [TJ . Initially 
the tagged particle is at the origin xt(0) = 0. Since 
particles do not pass each other, their order is clearly 
maintained, and the number of particles to the left and 
right of the tagged particle N is fixed. 

In Fig. [T] the straight line which starts at a; = is 
called the Jepsen line and follows x{t) = vt, where v is 
a test velocity [HI, |6| • We label the interacting particles 
according to their initial position, increasing to the right 
(see Fig. [T]). The tagged particle starts just to the right 
of the Jepsen line so at i = the label of the particle 
to the right of the Jepsen line is zero. In this system 
we have 2A^ -I- 1 particles. Hence initially we have N 
particles to the left of the Jepsen line and including the 
tagged particle iV -I- 1 particles to the right. 

Let a{t) be the label number of the first particle sit- 
uated to the right of the Jepsen line. According to our 
rules, at t = we have (5 = 0, and then the random 
variable a will increase or decrease in steps of -1-1 or —1 
according to: 

i) if a particle crosses the Jepsen line from left to right 
(5 — > (5 — 1 

ii) if a particle crosses the Jepsen line from right to left 
a — > ci -|- 1. Thus the counter a is performing a ran- 
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FIG. 2: In one dimension, the paths of a pair of particles 
in a hard core collision event can be represented by two non 
interacting particles which pass through each other, and then 
their label numbers are switched. 



dom walk decreasing or increasing its value +1 or —1 at 
random times. 

A collision between two hard core particles is repre- 
sented schematically in Fig. [21 In one dimension a hard 
core collision event is equivalent to two particles that 
pass through each other, i.e. non interacting particles, 
and then after the particles cross each other, the labels 
of the pair of particles are switched (see Fig. [2]). In- 
stead of relabeling particles after each collision, we let 
particles pass through each other, and then at time t we 
label our particles (or if we are interested only in tagged 
particle, locate the central particle). Operationally this 
means that in the time interval (0,t) we view the parti- 
cles as non interacting, and then find the particle with N 
particles to its right and N to its left, which is equivalent 
to the tagged particle in the interacting system. Hence 
the problem is related to the mathematical topic of order 
statistics [S^l as we will discuss briefly later. 

Following Jepsen Q and Levitt [H] we introduce the 
stochastic process for the non interacting particles a{€) 
where: 

i) if a particle crosses the Jepsen line from left to right 
a — > a — 1 

ii) if a particle crosses the Jepsen line from right to left 
a ^ a -f 1. The process a(t) is the same as the process 
a(i), in statistical sense. 

The event where a(<) = implies that the number 
of particles crossing the Jepsen line from right to left is 
the same as those crossing from left to right. Switching 
back to the interacting system, the event a(i) = means 
that the particle to the right of the Jepsen line at time 
t = is also there at time t. Hence the probability that 
ct(t) = is the same as the probability of finding the 
tagged interacting particle to be the first particle to the 
right of the Jepsen line. Therefore we will calculate the 



probability for a(t) = and with it obtain in the large 
N limit, the probability P(xt)Axt of finding the tagged 
particle in a small interval xtiXt + Axt- 

The random variable a is a sum of many random vari- 
ables 



N 



a 



(2) 



-N 



where 5aj is the number of times particle j crossed the 
Jepsen line from right to left minus the number of times it 
crossed the line from left to right. Clearly 5a j may attain 
the values —1 or if the particle started on the left, or 
1 or if it started on the right. Since we are interested 
in the large N limit we may neglect the contribution of 
(5ao , which is indicated by the prime in the sum [s^ . 

We calculate the probability Pat (a) of the random vari- 
able a. For that we designate Pll{xq'')^ the probability 
that particle — j, starting to the left of the Jepsen line 
at t = at Xq-* < 0, is found at time t to the left of 
the Jepsen line. Clearly for the corresponding trajectory 
5a^j ~ 0, since the particle crossed the Jepsen line an 
even number of times, or did not cross it at all. Similarly, 
Plr{xq'') is the probability to start to the left of the line 
and end to its right, and Prr{xq), Prl{xq) are defined 
similarly for particles starting on a;g > to the right of 
the line. 

In our non interacting picture, the motion of particles 
is independent, hence we can use random walk theory 
pol ] and Fourier series to find Pat (a). It is convenient to 
rewrite Eq. Q 



N 



(3) 



where Aaj — Saj 



Sa-j may attain the values 1, 0, —1. 
Each summand Aaj takes into account one particle start- 
ing to the left of the Jepsen line and one to the right. 

First consider = 1, that is one particle that starts 
on the left of the Jepsen line and one to its right. Then 
a = Aai — dai + Sa^i and as mentioned a = 1 or or 
— 1. Examples for possible trajectories are shown in Fig. 
[31 It is easy to see that 

( a^l Prl{xI)Pll{.x^^) 
Pn^i («) = < a = Pll{x^^)Prr{xI) + Plr{x^^)Prl{xI) 
[ a = -1 Prr{xI)Plr{x^^). 

(4) 

We define the structure function 



\{(t),Xf^- 



I Xq) — 



'^PLL{x^')PRL{xi) -f PLL{x^')PRR{xi) + PiflX (x^o) 



+ e-'f'PLR{x^')PRR{xi). (5) 



The structure function describes a single step in the ran- dom walk Eq. ^ in the following usual way: the co- 
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efficient of exp{i4>) is the probability that Aaj is equal 
one (i.e. PllPrl), the coefficient of exp{i(f>) = 1 (i.e. 
(j) — 0) yields the probability Aaj — 0, and similarly for 
exp(— i(/<). Since the summands in Eq. ([3]) are indepen- 
dent, Fourier analysis gives: 



PNia) 



1 

2^ 



(6) 




We average Eq. © with respect to the initial conditions 
xq which are assumed independent identically distributed 
random variable and we find 



{Pn (a)) - 
where from Eq. ([5]) 



1 

2^ 



d</)(A (0)) 



N -iad 



(7) 



FIG. 3: Trajectories crossing the Jepsen line. In the left 
panel one particle started on L and ended on L the other 
on R and ended on R. These two paths are assigned the 
probability Pll{xq^)Prr{xq) and yield in Eq. Q a = 0. 
Similarly the trajectories on the right panel correspond to 
Pll{xq^)Prl{xo) and they contribute a = 1. 



(A(<^)) 



-J^{Pll){Pb.l) + {Pll){Pbr) + {Plb){Prl) 



-e-'*{PLR){PRR 

(8) 



Here {Pij) is the probability of starting in i ~ L,R (rel- 
ative to the the Jepsen line) and ending in j — L, R and 
(• • •) denotes an average over initial condition. The prob- 
ability (Plr) is given in terms of the Green function of 
the non interacting particle g{x, xq, t) and the initial den- 
sity of particles situated initially to the left of the Jepsen 
line fLixo) 



where 

(Aa) = (Prl) - (Plr) 
\s expected, and the variance — ((Aa)^) 
'Prl){Prr) + {Plr){Pll)- 



2 



(12) 

(Aa)2 is 
(13) 



Using the central limit theorem we have the probability 
of zero total crossing of the Jepsen line, namely a = in 
the limit N oo 



PNia - 0) 



1 



\/27rA^(TAa 



exp 



N{Aa) 



(14) 



{PhRivt)) = _fL{xo) / g{x,xo,t)dxdxa 



(9) 



We see that for (Plr) we average over initial conditions 
in the domain (— L, 0) weighted by fhixo) (since the 
starting point is to the left of the Jepsen line) and also 
integrate over x in the domain {vt, L) with the weight 
g{x,xo,t) (since the end point is to the right of the line). 
Similarly 



/i?(a;o) / g{x,xo,t)dxdxo. (10) 



(PRRivt)) 



As mentioned fL{xo) [/_R(a;o)] is the PDF of initial po- 
sitions of particles that initially are at a:o < (a;o > 
0) respectively, hence fL{xo) — when a;o > and 
J^-j^ /L(a;o)dxo — 1, while fR{xo) is normalized and non 
zero only in {0,L). (Pll) a-nd (Prr) are defined simi- 
larly, and they satisfy (Pll) = 1 - {Plr), and (Prl) = 
1~{Prr)- 



III. DYNAMICS OF THE TAGGED PARTICLE 

We now apply the central limit theorem to analyze the 
sum Eq. ^ using Eq. Q when iV — >• oo. The small 4> 
expansion of the structure function Eq. ([S]) is 



This simple result, valid for a large class of Green func- 
tions and initial conditions, is suited for the investigation 
of a large number of single file problems. 

The probability PN{ct = 0) is according to our earlier 
discussion equal to the probability that the tagged par- 
ticle is in the interval vt < xx < vt + 6x, where 6x is the 
distance between the Jepsen line and the position of the 
second particle to the right of the Jepsen line. We as- 
sume that 6x is small relative to the typical length scale 
the tagged particle traveled. For example if the den- 
sity of interacting particles p is a constant, we assume 
1/p « \J {{xry). Hence we may replace vt xt to 
transform PN{ct ~ 0) to the PDF of the tagged particle 

p{xt) m. 

Making the replacement vt — )■ xt, using Eq. we 
have in the limit of large TV 



C exp < — ; 



PixT) ~ 
N[{PLRixT))-{PRLixT))] 



(A(0)) = l + z(Aa) 



1 



{{Aar)^' + Oi(b') (11) 



2 [{Pll{xt)){Plr{xt)) + {Prr{xt)){Prl{xt))] 

(15) 

where C is a normalization constant. In Eq. (fTSjl and 
what follows {Pij{xT)) is given by Eqs. (|9ll0p with 
vt Xt- Eq. ([15]) yields the PDF of the center par- 
ticle P{xt) in terms of {Pij{xT)) which according to Eq. 
([To)) depends on the free particle green function and the 
initial conditions. Thus the information contained in the 
non-interacting Green function is sufficient for the deter- 
mination of the single file diffusion of the tagged particle. 
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A. Thermal Equilibrium 

In the long time limit, and in the presence of a binding 
potential field, e.g. harmonic field or particles in a box, 
an equilibrium is reached. Then initial conditions do not 
play a role. For example {Plr{xt)) = P'i?' 



II (xt) with 



1 



cxp 



dx 



I XT \ 

where Z is the normalizing partition function 



Z = 



cxp 



Vjx) 

hT 



dx. 



Similarly 



PPixr) = ^ y _ cxp 



Vjx) 



dx. 



(16) 



(17) 



(18) 



In Eqs. (|16ll8p we used the steady state solution, 
liijit^ao g{x,xo,t) = exTp[—V{x)/hT]/Z, which is Boltz- 
mann's distribution suited for a system in thermal equi- 
librium. Using Eq. (|T5|) the position PDF of the tagged 
particle is liiRt^ao P[xt) = P'^^{xt) 



P'^'^ixT) - Cexp. 



N[{P'^^\xt))-{PI''{xt))] 
A{Pl\xT)){Pn'{xT)) 



(19) 

If the potential is symmetric V{x) — V{—x), e.g. parti- 
cles in a box or harmonic field, we have for not too large 
XT, PT'ixT) ^ 1/2, P^ixr) ^ 1/2 hence from Eq. ^ 

P'^'^ixT) ~ Cexp {-N [{P'^'^ixT)) - {PL{xT))f} . (20) 

Expanding the expression in the exponent in xt (since 
N is large) we find using Eqs. ()16ll8ll9p 



P^'^ixT) 



2VN 
.HiZ 



cxp 



47V 2- 



(21) 



where with out loss of generality we assigned V{x = 0) = 
0. Hence the standard deviation is 



^{xt) 



87V' 



(22) 



The same expression is found in Appendix A using the 
many body Boltzmann distribution, and integrating over 
all the particles except the tagged particle. 



iV particles on — e. This problem was solved already by 
Aslangul [26l |. and here we recover the known result us- 
ing our formulas. The Green function g{x, Xq, t) of a free 
particle is 



g{x,xo,t) 



exp 



(x-Xq) 

4£)t 



(23) 



where as mentioned D is the diffusion coefficient of the 
free particle. With the specified initial conditions we have 



{Prl {xt)) = lim / S{xQ-e) 



XT exp 



{x-Xq) 

ADt 



XT 



-dx. 



Similarly 



{Plr {xt)) = 1- 



V^TTDt 



dx. 



dxdxQ = 



(24) 



(25) 



When XT <C V^Dt we Taylor expand in xt to find 
{{Plr) - {Prl}? ^XTf/irDt, using Eq. m we re- 
cover the result in [291 



P{xt) 



1 

TT 



cxp 



N{xt) 



21 



TTDt 



hence 



{{xt)') 



TlDt 

~2N' 



(26) 



(27) 



The diffusion is normal, in the sense that the mean square 
displacement increases linearly in time. However the dif- 
fusion of the tagged particle is slowed down compared 
with a free particle, by a factor of 1/N which is due to 
the collisions with all other Brownian particles in the sys- 
tem. Clearly the approximation breaks down if one is in- 
terested in the tails of P{xt), since we used xt \/2Dt. 
Though clearly when N is large, the probability of find- 
ing such a particle is extremely small (i.e. use Eq. (|26|) 
P{xT = Vtot) - exp(-iV2/7r)). 



C. Formula for {{xTf) 



B. Simple Illustration 

We now consider the situation of particles free of a 
force F{x) — with open boundary conditions L — ?> cxd 
where initially all the particles are on the vicinity of the 
origin. More precisely, the tagged particle is initially sit- 
uated at Xt = 0, N particles to its right on e — 0+ and 



We now consider symmetric potential fields V{x) — 
V{—x), and symmetric initial conditions. The latter 
means that the density of particles at time t = to 
the left of the tagged particle, i.e. those residing in 
xq < 0, is the same as for those residing to the right, 
fR{xQ) = fL{~xo) (e.g. uniform initial conditions). In 
this case the subscript R and L is redundant and we 
use f{xo) = fR{xo) = fL{-Xo), where f{xo) = if 
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xq < 0, f{xo) > and /(xo)da;o = 1- From symme- 
try it is clear that the tagged particle is unbiased, namely 
{xt) — 0. Further, since N is large we may expand the 
expressions in the exponent in Eq. (jlSp in xt to obtain 
the leading term 

(Prl) - (Plr) = 



d 



[{Prl {xt)) - {Plr {xt))] \.^=oXt + 0[{xt)^ 



dxT 

Similarly for small xt 

{Pll){Plr) + {Prr){Prl) = 
2{Prr {xt)) [1 - {Prr {xt))] U^=o + 0{xt) 



(28) 
(29) 



To derive Eq. ((28)) we used the symmetry of the prob- 
lem, which means that the probability of crossing the 
point Xt = from left to right is the same as the proba- 
bility to cross from right to left {Plr)xt=o = {Prl)xt=o 
and similarly {Pll)xt=o = {Prr)xt=o by symmetry. We 
designate the reflection coefficient. 



n = {p, 



RR)xT=0 



(30) 



[or TZ = {Pll)xt=o] since it is the probability that a 
particle starting at xq > 0, is found in x > at time t 
when an average over initial conditions is made 



n 



L pL 

f{xo) / g{x,xo,t)dxdxo. 

JQ 



(31) 



A transmission coefficient is defined through T ~ 
{Prl)xt=o — {Plr)xt=o which is related to the reflec- 
tion coefficient in the usual way T = 1 — TZ. 

Turning our attention to Eq. (|28|) , from left-right sym- 
metry we have 



d 



{Prl {xt))\xt=o 



dxT 

Hence we define 



d 
dxj 



{PLRixT))\xr=0- (32) 



dxT 



{Prl (a;T))Ur=o 



where from Eq. (jTU]) 



f{xo)9{0,XQ,t)dxo. 



(33) 



(34) 



So r is the density of non interacting particles at a; = 
for an initial density f{xo)- Note that since {Prl{xt)) + 
{Prr{xt)) = 1 we have 



d 
dxi 



{Prr{xt))\xt=o- 



(35) 



Inserting Eq. ([55]) in Eq. 
have {Prl) - {Plr) - '2rxT - 



using Eq. (1321) we 
Inserting Eq. (j5n|) in 



Eq. (f^^ we find our main result: the probability density 
function of the position of the tagged central particle 



P{xt) 



1 



V2^((^t)2) 



exp 



{xt) 
'2{{xt)') 



where 



,, ,2, n{i-n) 



(36) 



(37) 



is the mean square displacement of the tagged particle. 
The single particle probability {Prr{xt)) gives TZ Eq. 
P0|) and r Eq. ()35|) which in turn yield the mean square 
displacement of the tagged particle Eq. ([57]) . We will 
soon use this equation to demonstrate a variety of phys- 
ical behaviors. However first we establish a connection 
between our work and the theory of order statistics. 



D. Order Statistics 

Generate n independent identically distributed ran- 
dom variables drawn from a PDF f{x) and arrange them 
in increasing order. Order statistics deals with the m-th 
observable m = 1, • • • , n among n observations taken in 
the increasing order, which is denoted Xm- The PDF of 
Xm, 4>{xm) depends on the PDF f{x), the sample size 
n and the order m. Let R{x) be the cumulative distri- 
bution of x, e.g. if the domain of x is —oo < x < oo, 
R{x) = f{x)dx as usual. Following well known re- 
sult [s^l, define x with 



R(x) 



m 



1 



then when n is large and m/n is of the order 1/2 



0(X„ 



const exp 



n{x — x)^ f'^{x) 



Hence the variance of : 



(a™)' = 



2R{x) l-R{x) 



R{x) l-R{x) 



(38) 



(39) 



nf^ (x) 



(40) 



which has some resemblance to Eq. (fS7|) . 

The problem of the motion of a tagged particle is math- 
ematically identical to the problem of order statistics in 
two cases: (i) in the presence of a binding potential and 
in the long time limit and (ii) when all the particles start 
on the same point. In both cases the single particle PDF 
g{x, Xq, t) of all the particles are identical since it is either 
independent of xq (case (i)) or we have a unique initial 
condition (case treated in sub section PlI B|) . For exam- 
ple in the presence of a binding field limf^oo g{x, xq, t) — 
e-xp[—V{x)lki,T]/Z which is independent of the initial 
position of the particle. Hence in equilibrium, to find 
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the center particle we may draw (2N + 1) random vari- 
ables from Boltzmann's distribution, and search for the 
center particle (which will give the position of the inter- 
acting tagged particle). Or, using the language of order 
statistics we have f{x) — exp[~V (x) / k^T] / Z . Using a 
symmetric potential V{x) = V{—x) and Eq. ([38|) we in- 
sert n = 2N + 1, and m ^ N + 1; hence, when N is large 
we have 



R{x) 



exp 



dx 



Thus, a; = 0. Then using Eq. (1^01) we find 



((^t)'> 



8N 



(41) 



(42) 



which is the same as Eq. (gj) (recall 1/(0) = 0). 

While Eq. (|3D| has superficially a structure similar to 
Eq. (|57|) they are different. Our TZ Eq. ([50)) is gener- 
ally not equal half, neither must it be close to that value 
[so Eq. ([38]) is not generally related to our problem]. 
In fact when t — > we must have {{xt)'^) — > which 
is found when limt_j.o7?. = 1 (see examples below). In 
the problem of motion of a tagged particle, the number 
of particles which can interact with the center particle 
is usually increasing with time (hence n is not fixed). 
For example consider the classical case of uniformly dis- 
tributed particles in infinite space and in the absence of 
forces. Roughly speaking particles at distances of the 
order of ^cff — VDt or shorter can influence the tagged 
particle motion via collisions. So in this case roughly 
-^cff = pV Dt particles participate in the process. In con- 
trast, if all particles are initially at the vicinity of the 
origin N particles participate i.e. influence the motion 
of the tagged particle. Interestingly this gives an argu- 
ment for the well known behavior of the mean square dis- 
placement of the tagged particle, with uniform density of 
particles, namely use Eq. (P7)) {{xtY) ^ ■nDtjlN and 
replace TV with iVog = psfDt to find ((xr)^) ^ \lDtlp 
which of course misses the correct numerical prefactor 
(see Eq. [50] below) . 



IV. PHYSICAL ILLUSTRATIONS 



In the limit L — cxd and iV — )• oo in such a way that the 
density p = N/L is fixed, we obtain single file diffusion 
in an infinite system, a case well studied long ago |3, § ■ 
The single particle Green function of a parti- 
cle in a box, with refiecting boundary conditions 
dg{x,Xi),i) / dx\^ ,j^ = is solved using an eigenfunc- 
tion expansion [38| 



g{x,xo,t) 



1 

2l 



n— 1 



2L 



L) 



2L 



(xo + L) 



exp 



With Eqs. ([3ll43l44|) we find 



E 



exp ( -D^t 



n2^2 ^ 

4L y 
(44) 



(45) 



«=i,Odd 



where the summation is over odd n. Kt t — Q, TZ = 1 
since all particles initially in (0, L) did not have time to 
move to the other side of the box, and limt_j.oo 72. = 1/2 
since in the long time limit there is equal probability for 
a non interacting particle to occupy each half of the box. 
Using Eqs. P4l43l44p we find 



1 

21 



(46) 



hence Eq. ([57]) gives the mean square displacement of 
the tagged particle 



n{t) [i-n{t)]L' 

N 



(47) 



The eigenvalues of the non-interacting particle determine 
the multi exponential type of decay of TZ(t) with time, 
which in turn determines the dynamics of the interacting 
tagged particle [43] . 

Let = Dir'^t/AL hence the limit 6^1 gives the 
short time dynamics before the particles interact with 
the walls, or equivalently the limit of an infinite system. 
The reflection coefficient is rewritten 



A. Particles in a Box 

Consider particles in a finite box extending from —L to 
L which was recently treated with the Bethe ansatz and 
numerical simulations by Lizana and Ambjornson [25j . 
The tagged particle initially at = has N particles 
to its right and N to its left. These particles are assumed 
uniformly distributed, hence 



/(^o) = 



< Xo < L 



otherwise. 



(43) 



1 



- 1 



1 



.i,Odd 



n=l,Odd 



(48) 

When (5 <C 1 we may replace the first summation with 
integration 



E ■ 

n=i.Odd 



1 



1 



dy = 



(49) 

where the factor 1/2 on the RHS comes from the 
summation over only odd n on the LHS. Using 
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FIG. 4: For particles in a box we show {{xtY) N / L of the 
tagged particle, in equilibrium versus A'^. Exact expression 
obtained in [S^ valid for all (dots) converges in the limit 
of A'^ — >■ oo to the behavior predicted by our theory (dashed 
Une) Eq. (|52|1 . For TV = 70 clear deviations between our 
asymptotic theory and the exact result in [2^ are found. 



Er=i,Odd V"^ = 7^2/8 we have 7^ 
and hence 



= 1 



VDt 



'Dt 



ITT p 



(50) 



This result was obtained in [1, [1], and it describes the 
dynamics of the tagged particle in a box, before particles 
have time to interact with the walls, namely when 5 is 
small. 

In the opposite limit of long times and finite systems 
we attain equilibrium, then limt_j.oo TZ = 1/2 and 



Since our formalism is valid only in the large N limit, 
comparison of our solution to an exact result like Eq. 
([55]) provides insight to the question of convergence. 

In Fig. m we plot limt^oo((a;T)^)-/V/L versus N us- 
ing Eq. ([551) comparing it to our Eq. (I52|) which gives 

limt^ao{{xT)'^)N/L — 1/2. Not surprisingly we see that 
the two results yield asymptotically the same result. The 
figure illustrates that even for N = 100 deviations be- 
tween exact results and the present theory are observ- 
able. 

In (25| numerical simulation of systems with N — 
1,10,70 (maximum of 141 particles) were performed, 
and favorably compared with the Bethe ansatz solution. 
These simulations were made for finite size hard core par- 
ticles whose diameter is A. To make comparison between 
our theory and simulation we scale the size of the box 
according to 2L ^- 2L - (2iV)A. In Fig. [5] the scaled 
mean square displacement versus t/Tf,q is shown, where 

2 

Teq = 4i /D. Some general features of our theory can 
now be discussed. First, at short times simulations show 
normal behavior {(xt)"^) = 2Dt, which is a trivial ef- 
fect: particles did not have time to collide and hence 
the tagged particle diffuses normally as if it is free. After 
particles start to collide, deviations from normal diffusion 
are observed (for N = 70 but as expected not for = 1). 
These are mainly due to collisions with other particles. 
Finally saturation due to the finite size of the system is 
found. While our theory shows the general trend of sim- 
ulations (say for = 70 and for not too short times) it 
is clearly not in perfect agreement. We argue that this 
is due to the small number of particles A^ — 70, since as 
we showed in Fig. |31 at least for large times and large 
A^ our theory gives the exact result. It would be nice 
to have simulations with point particles, with larger A^ 
and since there are three phases of the motion: normal 
diffusion, single file diffusion without the influence of the 
walls and finally saturation to equilibrium. Separation 
of time scales is needed to demonstrate these behaviors 
clearly. 



PixT) 



N 



hence 



lim ((xt)^> = L I2N. 



(51) 



(52) 



Eq. (ISTjl is a special case of the more general Eq. (PT|) 
since the single particle normalizing partition function 
for our example ys, Z = 2L. 

As mentioned, in [25| the Bethe ansatz was used to 
solve the problem of tagged particle motion in a box. 
Among other things an exact expression for the long time 
limit of ((xt)^), valid for all A^ was found [1^ 



lim {{xtY) = L 



N+l 



T{\)T[2 {N + 
T{N + 1)T{N ■ 



1) 



(53) 



B. Gaussian Packet 

Consider particles without external forces V{x) = 0, 
in an infinite system. Initially particles are spread with 
a Gaussian packet with width ^: 



f{xo) 



V2_ 



exp 



2e 



(54) 



for xo > 0. As before we consider the tagged particle 
motion, which is initially at xt = 0, with A^ particles to 
its left and A^ to its right. With the free particle Green 
function g{x,xo,t) Eq. (^5)1 we proceed to find {{xt)^). 
The reflection probability Eq. ((3T|) is 



dxr 



2 (°=o) 
e 



da; 



_ (x-xo)f 



^/4ZD^ 



(55) 
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FIG. 5: Motion of tagged particle in a box is shown for 
A'' = 1,10,70. Simulation results (squares) are taken from 
psi ] (see text). At short time simulations exhibit normal dif- 
fusion {{xt)'^) = 2Dt as illustrated by the straight dashed 
line which is a guide to the eye. Our theory is expected to 
work well in the large A'^ limit and when many collision events 
between tagged particle and surrounding Brownian particles 
took place. Hence agreement between theory (solid line) and 
simulation is reasonable at most only for A'' = 70 and not for 
too short times. 



10^ 
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10'- 



10"^ 



10" 



10^ 



10^ 



2Dt/^' 



Changing variables according to = {x — xo)^/4:Dt 
and using dimensionless parameter ^ = ^/y/2Dt, we find 



n = 



1 



dioe 



(x 
2« 



o)!Erf (io/\/2) 



(56) 



where Erf(io/V2)/2 = J^" e'y ^^dy/V^n is the error 
function |42]. Mathematica solves the integral in Eq. 
(l56l) and we find 



1 1 / V2Dt 

TZ = — \ — arccot 

2 TT \ C 



/2Dt 



(57) 



, namely 



For short times V Dt <^ we have TZ ^ 1 
most particles did not have time to cross the origin, while 
in the opposite limit limt^c>o ^ — 1/2 due to the sym- 
metry of initial conditions. The calculation of r Eq. ([M)) 
using Eqs. (|23l54p is straightforward 



(58) 



Inserting Eqs. (|57|58|1 in Eq. ^ we find 



1+ 



2Dt\ 



1 1 



r arccot 




This solution is shown in Fig. |5] with its two limiting 
behaviors 



/2Dt 



ttD 
2N 



N 
t 



short times 2Dt < 
long times 2Dt > 



(60) 



For short times the particles do not have time to disperse; 
hence, the motion of the tagged particle is slower than 
normal, increasing as t^^'^ which is similar to the single 
file diffusion with a uniform density Eq. (150)1 . Roughly 
speaking, for short times the tagged particle sees a uni- 
form density of particles with p = N/^. For long times 
we recover the behavior in Eq. ^7} since the scale of 
diffusion is much larger than ^. Hence if we start with 
a Gaussian or delta function packet we get in the long 
time limit similar behavior, as we showed. 



C. Particles in Harmonic Oscillator 

Consider particles in a harmonic potential V{x) — 
muP'x^ (2 where a; > is the harmonic frequency. The 
sing le particle undergoes an Ornstein Uhlenbeck process 
[38| and the corresponding single particle Green function 
is 



FIG. 6: Scaled mean square displacement of the tagged par- 
ticle with Gaussian initial conditions of the packet of particles 
exhibits a transition between short time {{xtY) oc t^^'^ law to 
{{xt)^) oc t behavior. Dashed lines are short and long time 
asymptotic behavior Eq. ((60)l. circles represent Eq. (|59)l . 



]{x,Xo,t) 



2ttDt (l - e-2*/^) 



: exp 



xae 



-t/r\ 



2Dt (1 - e-2*/r) 



, (61) 



10 



For short times t << t 



- -O— 3 -O— S -O— S -0— 9 -O— 9 -O— 9 -0- - - 



(67) 



10"^ 10'^ 10° 10^ 10^ 
t/T 



Such a behavior is expected, as well known for short times 
the diffusion process of the Ornstein Uhlenbeck process 
X ~ t^/^ is faster than the drift x ^ t and dominates 
the process, i.e. take i << t in Eq. and get the 

Green function of the force free particle Eq. Hence 
for short times we can neglect the external field, but use 
of course the Gaussian initial condition Eq. (j62|. Then 
the problem reduces to the short time behavior of the 
Gaussian packet, free of external forces (compare short 
time behavior in Eq. ([60)) with Eq. ((67)) when ^ — > ^th)- 
The long time limit of Eq. dm gives {{xxf) ~ tt£.IJ^N 
which is the thermal equilibrium behavior predicted more 
generally in Eq. (PT)) . Behavior of the scaled mean square 
displacement and its asymptotic behaviors is presented in 
Fig.El 



FIG. 7: Scaled mean square displacement of tagged particle 
in harmonic field Eq. (|66|) exhibits a transition between a 
short time {(xt)^) oc t^^^ law to saturation due to the binding 
field. The horizontal dashed line is the long time ((xt)^) = 
■K^fi^/AN behavior, the dashed line is the short time limit Eq. 
((67t , and the dotted dashed line is Eq. ((66)) . 



where (r) ^ — Dmcu^ / k^T is the inverse relaxation time. 
We assume thermal initial conditions 



= ,^ , ^ cxp 



'muj 

Using Eq. ((Ml) it is easy to show that 
_ 1 



(63) 



where the thermal length is ^th = V Dt = y/ muj'^ / ki,T . 
Note that one can write r — 1/Z which as we will soon 
prove is a general result valid for all potential fields satis- 
fying V{x) = V(—x), provided that the initial condition 
/{xq) is the thermal equilibrium . The reflection coeffi- 
cient [Eq. dni)] is 



n = 



;Erf(-^=|<i„, 



(64) 



where 77 = e^l'^ \J\ — e^^t/^. Using Mathematica 

- + -arccot f Ve2*/^ - l) . (65) 
1 - %/2i/(7ryr), and for 



7^ = 

For short times i/r << 1, 7?. - 
long time Tl^Xjl^ e-^/^ /tt. 

Using Eqs. p7l63l65p the mean square displacement of 
the tagged particle is 



1 



1 

r- arccot' 



(Ve2*/^ - 1) 



(66) 



D. Thermal Initial Conditions 

If we assume that initially the particles are in thermal 
equilibrium, our simple formulas simplify even more. If 
/(xo) = 2exp[-y(a;o)/fcfcr]/Z then r = l/Z. To see this 
notice that for symmetric potentials V{x) = V{—x) we 
have (?(0, — a:o; t) = 5(0, xq, t) and 



xo > 0. (62) yields 



exp 



exp 



khT 



g{0,xo,t)dxo (68) 



V{xo) 



z 



g{0,XQ,t)Axo, 



(69) 



where we assumed that the potential is binding so a sta- 
tionary solution of the Fokker-Planck equation is reached; 
i.e. the free particle is excluded. Therefore r [Eq. (l69l) ] 
is the probability of finding a non interacting particle 
at the origin, with thermal initial conditions. Since the 
thermal equilibrium density is the stationary solution of 
the Fokker Planck operator, r is time independent and 
equal to r = exp[— y(0)/fcbT]/Z. We can always take 
V{Q) = and then r = 1/Z. Examples for this behavior 
are the already analyzed cases of particles in a box with 
uniform initial density [Eq. ((46)) ] (since Z = 2L for that 
case) and similarly for particles in a harmonic potential 
with initial thermal density [Eqs. ((^ IS5)) ]. Hence using 
Eq. (127]) 



7^ (1 - 7^) z^ 

27V 



(70) 



For the symmetric and binding potentials under consid- 
eration, we have limj^oo TZ = 1/2 and hence Eq. ((70() in 
the long time limit yields the equilibrium behavior Eq. 
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E. Power Law Type of Initial Conditions 

Flomenbom and Taloni [23] considered particles free of 
external forces, hence the free particle Green function is 



g{x,xo,t) 



1 



exp 



{X - Xq) 

ADt 



(71) 



with initial conditions of power law type 

fixo)=B\xQ\-^ 0<xo<L (72) 

for < /3 < 1. Here the system size L — >■ oo and the 
reader should not confuse L with L. The normalization 
of the PDF of Eq. ^ yields B = {1 - (3)1^-^ In the 
limit /3 — )■ and L — )■ oo in such a way that p = N/L 
remains fixed, we anticipate the classical case of single 
file diffusion in the presence of a uniform density of par- 
ticles: (x^) oc t-^/^ [Eq. ([50)) ]. In the opposite limit of 
/3 — >■ 1 the particles are initially concentrated at the ori- 
gin, and we expect the behavior {x^) oc t [Eq. (P7|) ]. 
Thus < /3 < 1 bridges between these two known be- 
haviors, indeed one finds {{xt)^) c>c as shown in 
[13] ■ The latter is valid for times \/4:Dt ^ L as discussed 
below. This indicates that specially chosen initial condi- 
tion may control the qualitative behavior of the diffusion 
of the tagged particle. Here we analyze this case using 
our formalism finding analytical expressions for the mean 
square displacement. 

To obtain ((xt)^) all we need to do is to find r and TZ. 
Using Eqs. (I34l7ll72l) 



solving the integral 



(1-/3) 1 




y-^e~y'dy, (73) 



ADt / 



74) 



where r(a, z) is the incomplete Gamma function [42 
The following behaviors are found for short and long 
times 



(75) 




In the short time limit we took the upper bound in inte- 
gration in Eq. (f73l) to cx) so 

(1-/3) r -P ^y^. l-/3 r(l-/3/2) 
Cr ~ 1= — / y ' e ay — — 1= . 

VTT Jo yJ-K 2 

(76) 

Inserting Eqs. (|71l72p in Eq. (PT|) we find the reflection 
coefficient 

(77) 



In the limit of long times we get the expected behavior 
lim(_j.oo 7?, = 1/2 since then half of the particles are to the 
left of the origin and half to the right (in statistical sense) . 
Mathematica solves the integral in the last equation in 
terms of tabulated functions 

/;2/-^Erf(y)dy = 

[x [V^Erf (a:) + xE^/^ (x^)] - x^Y (l - f ) } , 

. 

where E„(z) = di cxp(— zi)/t" is the exponential in- 
tegral function. 

To analyze the short time behavior we use integration 
by parts in Eq. ((77]) and find 



7^ 



{^\ Erf {^\ ^ L^/^ yi-^e"^^dy 

\V4Di J \ViDi ) ^/tt Jo " » 

(79) 

Using asymptotic properties of the Erf function [421 ] , ne- 
glecting terms of the order of exp(— L^/4Z3t), Eq. ([75]) 
for y/lDt < L gives 



/4Dt 



P-1 



1-/3 



n - 1 



L 



(80) 



with 



dy 



r 1- 



2^ 



(81) 



Eqs. ()74l77l78p yield the mean square displacement 
Eq. ((371). The short time behavior y/ADt < L 



{(xt) 



1 



2iV(l-;3)2r (1-/3/2) 
while for long times y/ADt 3> L 



L 



1+/3 



{{xrf) 



nDt 

'2N' 



(82) 



(83) 



Thus after a long time the diffusion is normal ((x^)^) oc i, 
exactly like the case where all particles started initially 
on the origin Eq. ([?7|) . While for short times ((xt)^) oc 

i(l+/3)/2^ 

Taking the limit /3 — in Eq. ([5^ . we get the well 
known result of single file motion in uniform density of 
particles [1, [i] Eq. ((SU]), with the density p = N/L. In 
the opposite limit /3 — > 1 we have (x^) = for all 
times which is the same as in Eq. ([?f[) . Note that the 
limit /3 — > 1 must be treated with care, the short time 
limit and the /3 — >■ 1 limit do not commute due to the 
the 1/(1-/3)2 divergence in Eq. ([82]). When ^ ^ 1 all 
particles are centered on the origin hence we get the same 
behavior as in Eq. ([?7)) . 
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FERGUS RELATION: BEYOND BROWNIAN 
PARTICLES 



where the mean of the absolute value of the free particle 
motion is by its definition 



So far we considered the case where particles are diffus- 
ing according to the laws of Brownian motion in between 
collision events. What happens for other types of mo- 
tion? For example what happ ens when the underlying 
motion itself is anomalous 

Percus [l^ (see also 



4411) investigated a general 



relation between the diffusion of the tagged particle 
and motion in the absence of interactions (free motion) 
{{xt)"^) ^ {\x\)fi-cc/p- Such a relation was suggested 
for normal diffusion where we have (|a;|)free cx: t^/^, so 
{{xtY) oc t^/"^ and for a particle moving ballistically be- 
tween collision events hence (|a:|)tioe c< t and therefore 
when we turn on the interactions ((xt))^ oc t. This sim- 
ple relation between free particle motion and the mean 
square displacement of the tagged interacting particle, is 
expected to work for an infinite system, with a uniform 
density p of particles, and when external forces are zero 
[F{x) = 0]. Here we will derive the Percus relation from 
our formalism, for more general dynamics. 

Assume that the Green function of the non interacting 
particle can be written in the scaling form 



g{x,xo,t) 



1 



X — Xq 



(84) 



We assume G{y) = G{—y) and from normalization 
J'_^ G{y)dy = 1 . Here the free particle motion is anoma- 
lous, so that (|a;|)fi.cc oc for < 7 and not equal 
unity. For example the underlying motion might be sub- 
diffusive continuous time random walk (CTRW) [45-47] 
or fractional Brownian motion, where in the latter case 
G{y) is Gaussian and in the former G{y) is expressed 
in terms of Levy distributions (see some details below). 
We assume that moments of the process are finite. The 
constant has units of m^/sec'*'. 

To obtain the mean square displacement of the tagged 
interacting particle we calculate the reflection coefficient 
Eq. dSID 



7^ = = 
L 



dxo / G 
Jq 



Xo 



dx 



(85) 



Here we used uniform density of particle f{xo) = 1/L, 
and we will consider the limit L — > 00 with the den- 
sity p being kept fixed. Change of variables y = (x — 



Xo 



)lJK~^Vl'^ and using L I ^ K^V ^''^ — > c» we have 



7^ ' 



dxtx 







G{y)dy. 



(86) 



Integrating by parts, changing variables y = 
xqI ^JW^ril'^, and using the normalization condition, we 
have 



7^~ 1 



(1^1) 



free 



2L 



(87) 



(kl)frc 



\x\G 



dx 



Here we used the assumed symmetry G(y) = G{—y) 
which means that the underlying random walk is not bi- 
ased, as a result {\x\) = yG{y)dy. It is 
easy to see that r = 1/2L. In the limit where TV — > 00 
and L — T' 00 we find using Eqs. p7l87p the Percus for- 
mula for a general class of stochastic dynamics 



\'^\) free 

p 



(89) 



Clearly this equation gives a useful relationship between 
motion of a free particle and the same particle moving 
in single file when it is surrounded by identical particles 
whose density is p, 

As a simple example consider frac- 
tional Brownian motion where g{x,xo,t) — 
(^47rA%tT)-i exp [-{x - xq)^ / 4:K^t''] , where 
< 7 < 2. Then Eq. ^ gives 



P\^ 



(90) 



Hence if the underlying motion is ballistic 7—^2 the 
motion of the tagged particle undergoing single file dy- 
namics is normal with respect to time ((cct)^) oc t. For 
a CTRW particle in the continuum approximation the 
non-interacting single particle green function is governed 
by the fractional diffusion equation [47l - l49j 



—g{x,xo,t) = Kj—g{x,xo,t) (91) 

with < 7 < 1. From this equation it is easy to find (HO] 
(|a;|)froc = Jk-,V/T{1 + 7/2). Hence using Eq. 



pF (1 + 7/2) 



(92) 



We see that independent of the mechanism of the under- 
lying anomalous diffusion (i.e. fractional Brownian mo- 
tion, or CTRW) we get for the tagged particle ((xr)^) ^ 
V/"^. Further our general results show that the Green 
function of the tagged particle is Gaussian even though 
the Green function of non-interacting CTRW particles is 
highly non-Gaussian. 

Warning: One should be careful when applying our 
results to the CTRW model, namely Eq. ([^ should not 
be abused. One mechanism of anomalous diffusion are 
the power law waiting times of Scher and Montroll [i^ 
which yield slow dynamics, as captured by the CTRW 
[ii,!!^ and the fractional Eq. (EU "SC!]. Single file diffu- 
sion with such sub-diffusive motion as the starting point 
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was considered recently theoretically and with simula- 
tions in [13, HI] . It should be noted that the meaning 
of collision in such a model should be taken with care. 
For a random walk on a lattice with waiting times on 
each lattice point, one can envision several collision rules. 
For example one might allow two particles to occupy the 
same site at a given time, or one may consider a mech- 
anism were a particle once hopping into a trap already 
occupied will eject the particle previously residing in the 
trap. Or a particle is allowed to jump only into an empty 
site, as is usually assumed. These type of collision rules 
might yield behaviors different than ours. For example a 
particle stuck with a very large sojourn time, might be 
ejected by another particle, hence one can imagine a situ- 
ation where some form of interaction causes the particles 
to move faster. In our work the collision implies that we 
can let particles go past one another as if they were non 
interacting (so on a lattice two particles may occupy the 
same point at the same time) and eventually we look for 
the center particle. Even more interesting will be to in- 
vestigate interacting particles in systems with quenched 
disorder, since the latter, when strong enou gh, is known 
to lead to non Gaussian sub-diffusion [isl. |46|. There 
the simple formula Eq. is generally not expected to 
hold. Indeed in [s^] single file motion of a tagged particle 
in the Sinai model was considered, the results are much 
richer when compared with Eq. ([M)) . 



problem onto a non interacting problem using the Jepsen 
line. The motion of the tagged particle belongs to the 
general problem of order statistics. The problem reduces 
to considering a list of 2N -\- 1 random variables and find- 
ing the distribution of the variable which has N variables 
smaller than it and N larger corresponding to center par- 
ticle (note that the right most particle we will have an 
extreme value problem). Classical theory of order statis- 
tics deals however with the case where all the random 
variables have identical distributions. In contrast in the 
exclusion process under consideration, particles have non 
identical distribution. Thus except for two cases: i) all 
the particles initially on the same position and ii) equi- 
librium state, the problem deals with non identically dis- 
tributed random variables (since the initial condition are 
non identical). While we treated the problem of sym- 
metric potential and symmetric initial condition for the 
center particle in detail, it is left for future work to con- 
sider non-symmetric potential fields, non symmetric ini- 
tial conditions, and the dynamics of the particle in the 
tails of the packet. We believe that the methods devel- 
oped here with some modifications can treat these cases 
too. 

Acknowledgment EB is supported by the Israel Sci- 
ence Foundation, and RS by the NSF. We thank Ludvig 
Lizana and Tobias Ambjornson for sharing with us their 
data in Fig. [5] and for providing helpful insight. We also 
thank M. Lomholt and A. Taloni for useful discussions. 



VI. DISCUSSION 



The one dimensional problem of motion of a tagged 
particle interacting via hard core interactions with other 
particles was solved using the Jepsen line. The formalism 
we developed treats both Brownian and non Brownian 
motion in between collision events, is suited for rather 
general external fields acting on the particles, for open 
and closed system, and handles also different types of 
initial conditions. Following others we have mapped the 



VII. APPENDIX A 

In this Appendix we use Boltzmann's distribution for 
the interacting system to find the PDF of the tagged 
center particle in equilibrium. The multi dimensional 
PDF for 2N + 1 interacting particles, in the presence 
of an external binding field V{x), with V{x) = V{—x), 
acting on all of them is 



P(a;_Ar, • • • ,a;_i,a;o,a;i, • • • ,a;Ar) = exp 

(93) 

where Z2N+1 is a normalizing factor, and 0{x) is a step function: 0(x) = if x < 0, 0{x) — 1 for x > Q. The center 
tagged particle is x^ = xt- To find the PDF of xt in equilibrium, which we call P'"^{xt), we must integrate Eq. ([M)) 
over all coordinates besides xo — >■ xt 



N 

E 

3 = -N 



hT 



{X^N+I - X-n) S ix_N+2 - X-N+l) ■■■0 (xn - ^JV-l) 



kuT 



3 dx-AT J^l^ dx-N+i ■ ■ ■ j^'^^ da;_i exp 



f°° dxi f°° dx2 ■ ■ ■ r°° dxN exp 



(94) 



kbT 



We rearrange the integration limits, as explained in Fig. [5] 



XT PXT 

dx^N / dx- 



N+1 ■ ■ ■ 



XT PX-N 

dx-N I dx- 



■N+l ■ ■ ■ ■ 



Hence we can use 



XT pXT ^ 

dx-M / dx- 



-N+l • • • — 



XT rXT rXT rx-N 

dx^N I dx-N+1 ■ ■ ■ + / da;_Ar / dx_ 



Repeating this procedure we rewrite Eq. (IMl) as 



ET exp 
P'^'i {xt) = Nor <( / dx 



N+l 
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(95) 



^ rXT rXT 
- I dxN I dx 




dx- 



exp 



V(x) 



N 



exp 



V(xt) 



z 



-co 

(96) 



(97) 




FIG. 8: Left panel integration in the domain — oo < X-n < 
XT, x-N < x-M+i < XT is equivalent to integration — oo < 
X-M < Xt, —oo < X-N+i < x-N (right panel). 



where Nor is a normalization constant and Z is defined 
inEq. (HID. 

Thus Eq. describes a problem of order statistics 
which is extensively investigated by mathematicians, as 
mentioned in Sec. HTlDl P'^'^ixT) Eq. ^ is the PDF of 
the random variable which has exactly N random vari- 
able larger than it and N smaller. In this sense we have 
transformed the problem to a non-interacting system, 
similar to the non interacting picture in the main text. 
The information contained in the single non-interacting 
particle, i.e. the single particle Boltzmann distribution 
exp[—V{x)/kbT]/Z is all what is needed for the calcula- 
tion of the position of the tagged particle. 

Using the symmetry V{x) = V{—x) we have 



-da; — 



1 



-da;, 



(98) 



and 



-da;. 



J X T ^ ^ JO 

Using Eqs. (|98I99|) . we rewrite Eq. (|92l 



P^'i (xt) = Nor 2 

f r /r^e- 

exp < In 



1-4 ^ ^ da; 



In 

fXT 



the limit N 



(99) 



(100) 



with 



oo only Xt 
jQ'^exp[—V{x)/ki,T]dx/Z <C 1 will have a measur- 
able contribution to P^'^{xt), since if this condition is 
not satisfied, the value of P'''^{xt) is exponentially small 
in N. Expanding the In in Eq. (llOOp we have 



P°i(a;T) oc exp 



hT 



exp 



Z 



(101) 

Since N >> 1 we expand J^, ^ e:iip{—V{x)/ki,T)dx = xt 
where we used V{Q) = 0. We use V{xT)/kbT <C 
N{xtY /Z"^ which holds in the center part of the PDF 
of the tagged particle (since >> 1, V'(O) = 0, and 

where C 



4N{xt) 

Z2 



V{x) is analytic) P''^{xt) ^ Cexp 

is a normalization constant. This final approximate re- 
sult is the same as Eq. (PT|) . justifying the tricks used to 
derive our main results. 
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